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Abstract 

The Wiener-Hopf factorization is obtained in closed form for a phase 
type approximation to the CGMY Levy process. This allows, for the ap- 
proximation, exact computation of first passage times to barrier levels via 
Laplace transform inversion. Calibration of the CGMY model to market 
option prices defines the risk neutral process for which we infer the first 
passage times of stock prices to 30% of the price level at contract ini- 
tiation. These distributions are then used in pricing 50% recovery rate 
equity default swap (EDS) contracts and the resulting prices are compared 
with the prices of credit default swaps (CDS). An illustrative analysis is 
presented for these contracts on Ford and GM. 
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1 Introduction 



The equity default swap contract has recently been introduced in the financial 
markets. The contract pays a recovery proportion of a notional amount at the 
time of the occurence of a specific equity event, up to a maturity, that is typically 
5 years. The equity event is defined as the first time the stock price drops below 
30% of the price prevailing at contract initiation. The distribution of this first 
passage time is then of critical importance in the valuation of the contract. 
These distributions are known for the geometric Brownian motion model and 
results have also recently been obtained (Davydov and Linetsky (2001), Campi 
and Sbuelz (2005)) for the constant elasticity of variance (CEV) model. The 
prices of European options under these models, however, are not consistent 
with observed market option prices and such observations call into question the 
resulting swap prices. 

Observed market option prices arc more closely matched by a variety of 
pure jump Levy process models for the evolution of the logarithm of the stock 
price. Here we shall focus attention on the four parameter CGMY Levy model 
introduced by Carr, Geman, Madan and Yor (2002). First passage times for 
such Levy processes require knowledge of the distribution of the supremum 
and infimum of the process over a fixed time interval. The product of the 
Laplace transforms of these distributions is known by the famed Wiener-Hopf 
factorization. Identification of the law of the infimum and supremum then 
flounders on the inability to analytically perform the factorization. Again, it 
is known that if the Levy process has only one-sided jumps, say downwards, 
then the law of the supremum is known and one may then use the Wiener-Hopf 
identity to deduce the law of the infimum. Such strategies have been effectively 
pursued in Rogers (2000), Novikov, Melchers, Shinjikashvili and Kordzakhia 
(2003), Khanna and Madan (2004) and Chiu and Yin (2005). Processes with 
one sided jumps are dominated by those with two sided jumps in their ability 
to explain option prices. 

For Levy processes with two sided jumps the double Laplace transform in 
time and the level of the infimum and supremum have been derived by Nguyen 
and Yor (2002) for some Levy processes using bivariate integral representations. 
However, the numerical inversion of these transforms involves four integrals and 
is computationally quite involved. Here we consider an alternative strategy for 
obtaining the first passage distribution of two sided jump Levy processes. To 
this end, we use phase-type distributions, defined as the time until absorption 
of a finite state continuous-time Markov process (See Asmussen (1992)). It is a 
well known fact that this class comes close to being complete for obtaining first 
passage probabilities of the type we consider. 

We employ results from Asmussen, Avram and Pistorius (2004) to obtain 
closed forms for the Laplace transform of first passage distributions for Levy 
processes when the Levy measure has a phase type distribution. In addition to 
providing tractable computational schemes, phase- type distributions have the 
advantage of forming a dense class so that by increasing the size of the state 
space, one can in principle get arbitrarily close to a target density and thereby 
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get arbitrarily good approximations of the first passage probabilities. 

A particularly important subclass of phase-type distributions is formed by 
the hyperexponential distributions, i.e. finite mixtures of exponentials. The 
underlying Markov process then chooses some initial state, stays there an expo- 
nential time with rate depending on the state and goes directly to the absorbing 
state. This subclass is appealing for the CGMY Levy case as the target density 
is completely monotonic, i.e. a possibly continuous mixture of exponentials, 
so that by a discretization procedure it is a limit of (possibly scaled) hyper- 
exponential densities. In practice, the fitting can be done either by maximum 
likelihood (Asmussen, Nerman and Olsson (1996)) as then one wants to capture 
the shape of the target density, or by minimizing a distance measure paying 
specific attention to specific features of interest. In particular, in first passage 
problems a good fit in the tail is crucial, and hyperexponential approximations 
to Pareto tails have been obtained in this way (Feldman and Whitt (1998)). 
This is the path we follow in this paper. We approximate the Levy density of 
the CGMY process by a hyperexponential distribution and we then use the first 
passage time of the approximating process to determine the Laplace transform 
of the first passage time distribution for the approximating Levy process. One 
dimensional Laplace transform inversions using the methods of Abate and Whitt 
(1995) then give us both the survival probability and the first passage density 
for this approximation. One may wish to evaluate the quality of this approx- 
imation to an exact result for the CGAIY process. The latter is however not 
computationally available at this time. One may obtain alternative approxima- 
tions by either solving partial integro differential equations or by Monte Carlo 
simiilation. The simulation would require the development of appropriate im- 
portance sampling methodologies. We therefore leave these questions for future 
research. The same considerations apply with respect to constructing a priori 
bounds on the error with respect to the exact result. 

From the first passage density and survival probability one may determine 
the quote on an equity default swap contract. We estimate the CGMY pa- 
rameters from the prices of options on Ford and GM over a three year period 
and extract daily estimates of quotes on equity default swap contracts using the 
above hyperexponential approximation. We compare these prices with market 
quotes on the CDS rates and we observe that these two separate sets of prices 
for a credit event, indeed, correlate well. 

The outline of the rest of the paper is as follows. Section 2 presents the 
details for the construction of the Laplace transform of the first passage time 
for an approximation to the CGMY Levy process. Section 3 describes how 
this transform is used in quoting on equity default swap contracts. Section 4 
presents the results of daily calibration of the CGMY model to market option 
prices. Section 5 compares the resulting equity default swap prices with their 
credit default swap counterparts. Section 6 concludes. 
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2 The CGMY First Passage Time 



We wish to determine the first passage time of the CGMY Levy process to 
various levels using a phase type distribution to approximate the process. For 
the Laplace transform of the first passage time of the approximating process, 
i.e. a process with a phase type Levy measure, we follow Asmussen, Avram 
and Pistorius (2004). Given that the Levy measure is completely monotone we 
may employ a very special phase type distribution that approximates the Levy 
measure as a mixture of exponentials. The phase type distribution we employ 
therefore has a simple form. Starting from an initial state we have an entrance 
into a number of states where we stay until absorption. There are no transitions 
between states. Suppose the Levy measure has the form 



i=l 



+ ^6,e-^^l^ll.<o. (1) 



(2) 



Our phase type distribution for the positive jumps has n + 1 states with the 
generator matrix 

/ -ai \ 

— Q!2 



with unnormalized initial state probabilities proportional to ^ for the states 
i = 1, - ■ ■ ,n and the vector of absorption rates into the final state is 

t = {ai,a2, ■■■ , an). (3) 

The characteristic exponent for the approximation now takes the special 
form 

k(s) = logi^ [exp (sXi)] 



+ s J 2 



n m 1 

= E^;a. = e|^ (4) 

_i_ ai _ bj 
^ 1. ' " TIT- 

Where the value of ct^ will depend on the diffusion approximation we incorporate 
for the small jumps. This is described in greater detail later. The drift /x is 
chosen to make the overall drift on the stock price be the interest rate less the 
dividend yield. 
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Let Tx be the first passage time over a level a; > 0, the Laplace transform of 
Tx at the transform argument a 

E [e"""^"] 

is obtained in terms of the roots, with positive real part, of the equation 

k{s) = a. (5) 

We observe that one may write 

as a ratio of two polynomials and so wc seek the roots of the polynomial equation 

p{s)=aq{s). (7) 

The polynomial q{s) is 

n m 

and 

+ MS - (A+ + A_) - A+ V TT+a,^— + A. V njPj -— 

(9) 

We observe that k(0) = and the poles of k(.s) arc exactly equal to the 
eigenvalues of — T'^+^ and those of T^"^ (and also the roots of q{s) = 0).We 
note that p — ag is of degree n + m + 2if(T>0 and of degree n + m + 1 if 
(7 = 0, /i > 0. In the case /x > 0, tr = 0, n{s) — a has n+1 distinct positive roots 
and m distinct negative roots. If cr > k(s) = a has n + 1 distinct positive 
roots and m + 1 distinct negative roots. Finally for ^ = 0, ct = there are n 
positive and m negative distinct roots. 

Let fc+ be the number of positive roots. Suppose the positive roots are p^, 
i = 1, - ■ ■ ,k+ then letting 

Mt = sup {Xs V 0) (10) 

s<t 

and writing 

P (M^^a) €dx) = / oe-"*P(Mt e dx) (11) 

(e(a) is an independent exponential random variable with mean a~^), we have 
that (for s with n{s) < 0) 
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/•OO 

Jo 

poo 

= P (Me(„) = 0) + / e^-P (Me(„) G rfa;) 

/•OO 

= 0+(-oo)+/ e^^P (M,(„) G dx) . (12) 
This function may be expressed in terms of the positive roots Pj by 

= det(-s7-r) n'i(-pJ 



(13) 



where p^, ■ • -p^.^ arc the positive roots of k{s) = a and where T = diag(—ai) 
is the restriction of r+ to the nonabsorbing states. If cr = the number of 
positive roots k+ is equal to the number of positive phases n, k+ = n and we 
note that 

k+ 

<t>t{-oo) = l[^ (14) 

and < (^+(-oo) < 1. 

If cr > then the number of positive roots is A;+ = n + 1 one more than the 
number of positive phases n : Pi, - ■ ■ Pn+i and we have the formula 

furthermore, now (/>^(— oo) = lims^_oo 4'a{~^) zero. 

Hence in the case that there is a Brownian component present (cr > 0) 

0+(-«3)=O. (16) 

By performing a partial fraction decomposition we can perform the Laplace 
inversion in s analytically: Note that we can write 



4,+ {s)-^ti-oo) = J2^t^^ (17) 

i=l ^ Pi 

where 

At = ^-^^-''^^^ . (18) 
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The Laplace transform can be inverted explicitly to find for a; > 

k+ 

P (M,(„) <x)=P (Me(„) = 0) + - e-".-). (19) 

It follows from the fact that 

P(M,(,) < oo) = 1 (20) 

that 

1-^A+=P (M,(„) = 0) = 0+(-oo). (21) 

i=l 

In particular we get that 



P(Me(a)>a:) = l-P(Me(a) =0) -^^+(1-6-"-^) 

i=l 

= l-EA+-^(^e(a)=0)+^A+e-''.- 

j=l i=l 

fc+ 

= YlA+e-P'''. (22) 



Now since 



i=l 



P (Me(,) > x) 



f-OO 

/ ae-"*P (Mt > x) dt 
Jo 

POO 

= a e-"*P (T^ < t) rft (23) 
Jo 

we find the Laplace transform of P (T^^ < t) 

/ e-"*P (T, < i) rfi = - V A+e"''*^ 
■^0 i=l 



(24) 



where = Pi{a) and = [a) depend on a. 

We employ the methods of Abate and Whitt (1995) to obtain first passage 
probabilities and the first passage density. We next consider the details for 
the diffusion approximation and the phase type approximation of the CGMY 
process. 

2.1 CGMY details 

The CGMY process was introduced in (Carr, Geman, Madan and Yor (2002), 
see also Koponen (1995), Boyarchenko and Levendorskii (1999, 2000)) and is a 
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pure jump Levy process with the Levy density 



M^)=^^l->o + ^^wl-<o- (25) 

For this completely monotone Levy density of the CGMY process we recognize 
that 

1 r<x> „,Y „—ux 

du (26) 



lo ni + Y) 
and consider the approximation scheme 

JV — 1 Y —ux/ \ 
1 _ Ul e "'^(Mt+l - U^) 



T{l + Y) 



(27) 



For a prespecified sequence of exponential decay coefficients ui that cor- 
respond to reasonable levels of mean jump sizes under the single exponential 
model. The specific values for Ui were obtained using a least squares algo- 
rithm that minimizes the sum of squared errors between the left hand side 
of ([77)1 and the right hand side of l|?7|) evaluated at the x points starting at 
X = .25 and increasing to a; = 5 in steps of .025. As starting values we employed 
.5, 2, 5, 10, 20, 40 and 100. Our application works with the value of Y fixed at 
.5 and hence the approximation of equation ([27l) developed and subsequently 
used is independent of all parameter variations. The resulting values for Ui were 
.1940, .5982, .8434, 1.1399, 1.5308, 2.1211, and 3.4055. 

We then approximate the CGMY Levy density on the two sides by 

N-l 

k+{x) = ^ c,e-(*^+"-)- (28) 

i=l 

^-ix) = ^d.e-(«+"')^ (29) 

- d.- r(i + y) ("^"^ 

We then have in terms of our general discussion above 

M + m; p.^^G + m (31) 

Alternative but somewhat related procedures for fitting hyperexponentials 
to general densities (rather than Levy measures as here) have been considered 
in Feldman and Whitt (1998), Andersen and Nielsen (1998) and Asmussen, 
Jobmann and Schwefel (2002). 
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2.2 First Passage to a low level 

For the first passage to a low level we may work with the first passage of —X{t) 
to a high level. In this case we reverse the roles of a, /3 and write 

a, = G + Ui (32) 

/?, = M + Ui. 

2.3 The small jump diffusion approximation 

Above we discussed a procedure to approximate a CGAdY process X by a Levy 
process Xi with density ([1]) and to improve the approximation we would like 
now to approximate the difference by a Brownian motion. 

A way to refine the above approximation of X by Xi is to approximate the 
process of small jumps by a Brownian motion. This process of small jumps 
denoted Zf is obtained using the Levy density of X restricted to (— e,e). The 
process Zf is approximated by a Brownian motion with variance 

= a^ie) = [ x^v{dx). (33) 



In Asmussen and Rosinski (2001) it is shown that, if v has no atoms, Zf 
weakly converges to a Brownian motion if and only if 

1 

lim — / X v(dx) — oo. (34) 
For a Variance Gamma process this diffusion approximation fails as 

I r _ 

lim - xe "■^dx < oo. (35) 
For a CGMY process with Y > Q this diffusion approximation is valid, since 



limi / x^-^e-'^^'dx^oo (36) 

(as v{dx)/dx > const/\x\^~^^ in a neighbourhood of the origin, see Example 
2.3 in Asmussen and Rosinski (2001)). 

In our CGMY process approximation we truncated the GGMY Levy density 
at ui — .5 on the one hand but, on the other hand, the approximating densities 
k^{x), k-{x) are taken to start at a; = 0. Therefore we apply the diffusion 
approximation to 

k{x) = (e-^^7xi+^ - k+{x)) lo<.<„, + (e-«l-l/|a:|i+^ - fc.(a;)) l_„,<,<o. 

(37) 
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3 Equity Default Swap Pricing 



We essentially follow the logic for the pricing of credit default swaps, replac- 
ing the required survival probabilities and default time densities with the first 
passage complementary distributions F„(C, G, M, Y) and the first passage time 
density fn{C, G, M.Y) as computed by Laplace transform inversion using the 
approximation methods described above in section 2 for n days or time 360t = 
n. 

Like the credit default swap, the equity default swap contract is viewed in 
two parts, one describing the receipt side of the cash flows and th(^ other the 
payment side. There is a notional amount M associated with tlie receipt side 
and the actual level of cash flows received is M times one minus the recovery rate 
R on the occurence of the equity event. Our calculations employ a recovery rate 
of 50%. These funds are received on the time r of first passage of the equity to 
a level below the specified barrier. We employ two barriers in our calculations, 
these are 50% and 30% of the price at initiation. 

Against this stream of receipts the equity default swap holder makes periodic 
coupon payments, typically until the equity event or the maturity whichever is 
comes first. At the first passage event date we subtract from the receipts the 
accrued coupons from the last coupon date before r to this date. We denote 
time measured in days by n and in years by t. For each day n we define the 
function ^(n) that gives the number of days between the end of day n and the 
last day on which a coupon was paid. We also denote by N the maturity of 
the contract in days while T is the maturity in years. By convention we take 
{N,n) = 360 (r,t). 

Let A(n) be the first passage event indicator function that takes the value 
1 if the first passage of equity has occurred on or before day n and is zero 
otherwise. Further, let kx denote the annual coupon rate or the equity default 
swap rate quoted on the contract for maturity T.The cash flow receipts on day 
n, Rn are then written as 

R. = (A(n) - A(n - 1)) (m(1 - R) - M^lh__m'^ . (38) 

Let NP denote the number of coupon payment dates and let npj be the day 
number of the the j*'* payment date. The j*'' payment Pj occurring on day npj 
is then given by 

P. = ^^^^^^l5^^(l-A(np,)), (39) 

where the multiplication by the complementary flrst passage event indicator 
function recognizes that no coupon payments are made after the first passage 
event. 

For the valuation of claims we employ a discount function B„ that gives the 
present value of a dollar promised on day n. There are a variety of approaches 
to constructing such discount functions. 
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The random present value of cash flows to the equity default swap contract, 
V{T) is then given by 

N NP 

y (T) = ^ i?„B„ - ^ P, B„p^ . (40) 

n=l j=l 

We are supposing here that conditional on the outcome of the first passage 
event the expected values of future dollars are the same. This is equivalent 
to supposing that the first passage events of single names in the economy do 
not contain any information about macro movements in interest rates or that 
interest rate evolutions are independent of the first passage process. We leave 
for future research the modeling of joint evolutions of interest rates and first 
passage times. 

The equity default swap quotes in markets are set at levels consistent with 
a zero price at the initiation of the swap contract. Hence we have that under 
the risk neutral measure 

£:0 [V{T)] = 0. (41) 
From our risk neutral distribution for the first passage time we employ 

E'i [(1 - A(n))] = G, Af, Y). (42) 

Def 

The probability of first passage on a particular day n may be approximated by 
the density of default times the length of the day in years. 

eQ [(A(n) - A{n - 1))] - ^„(C, G, M, Y) 

Def 

« f„{C,G,M,Y)^. (43) 
We may then compute the value of the equity default swap quote as 

N 



eQ[V{t)] = J2 



n=l 
NP 



Mil - R) - - 



B„7r„(G,G,M,y) 



360 

E ^^fe^(-^^^-^.-^) ^„^^T-„,^ (c, G, M, Y). (44) 

Setting this value to zero and solving for the kr we obtain the CGMY equity 
default swap pricing model 

MG,G,M,Y)- il-R)j:t^B.n.iC,G,M,Y) 



Ef/i ^""''.^o'-'^ B^P^Fnp, (G, G, Af, Y) + ELi ^B^n^{C, G, Y) ' 

(45) 

Equation (j45|) provides us with a four parameter equity default swap pricing 
model and the parameters may be estimated from market prices of options on 
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the underlying name. In this way we obtain option impHed quotes for the equity 
default swap rates. 

By way of a stylized example we set the interest rate yield curve at a flat 5% 

level and used the values C = .5, G = 2, M = 10, and Y = .5 with recovery at 
50% and obtained the EDS prices for maturities of 1, 3, and 5 years at 161.97, 
336.65 and 439.54. 

4 CGMY Calibration to Option Data 

We obtained daily data on 5 year CDS rates on Ford and GM over the period 
25 February 2002 to 25 February 2005. There wore 696 trading days in this time 
interval. The credit default swap rates on these names saw a sharp increase over 
this period before they finally came down again. The mean CDS rates over this 
period were 308 for Ford and 249 for GM. The standard deviations were 132 and 
78 respectively while the minimum rates were 155 and 115 respectively. The 
corresponding maximal rates were 720 and 480. 

We employ the above algorithm to compute option implied equity default 
swap rates with barriers for the equity event set at the typical value of 30% 
of initiation price and a five year maturity. For this purpose we calibrate the 
CGMY model to European option prices for each day on each name. Equity 
option prices on Ford and GM traded on the New York Stock Exchange are 
typically available for American options. We mitigate the American feature by 
first employing out of the money options. Second we infer implied volatilities 
from American option prices using the Black Scholes model and then construct 
the corresponding European prices by the Black Scholes formula. The data is 
obtained from OptionMetrics and is available at WRDS the Wharton Research 
Data Service. 

We note that the case F = is the variance gamma model, that successfully 
calibrates any of the maturities. We are also aware that we may calibrate equally 
well with any specific value of the parameter Y. For our approximations to have 
a diffusion component we require F > and the structure of the approximation 
is determined for any fixed value of Y. We therefore froze the value of the Y 
parameter at Y = .5 and calibrated the parameters CGM for this frozen value 
of Y. The calibrations were done using maturities between 1 and 2 years, by 
minimizing the root mean square error between market and model European 
option prices. We report in Table 1 summary statistics of the CGM parameters 
for both the companies. 



TABLE 1 





FORD 






GM 








C 


G 


M 


C 


G 


M 


Median 


.6506 


1.9458 


11.0187 


.2171 


1.0084 


5.8031 


25% 


.3661 


1.3066 


9.9309 


.1664 


.6690 


4.7486 


75% 


1.0895 


4.0969 


11.3522 


.5582 


2.7802 


11.5872 
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Figure 1: Ford CDS in Blue, EDS in Red 



5 Equity Default Swap Rates and the CDS Rates 



We employed the calibrated CGM model with y = .5 to determine the equity 
default swap rates for the typical barrier of 30% for the payout on the equity 
event. We present graphs of the implied EDS and CDS rates for both companies 
over the 2002-2005 period evaluated every five days. The data on CDS rates 
for the five year maturity are readily available from Bloomberg. 

The results of regressing the CDS rates on the EDS rates are as follows. We 
observe that there is a strong link between between the prices obtained from 
these separate markets. 



FORD 
Const 
value 99.1820 
t — value 13.62 



0.2913 
33.37 



Slope 



i?2 

0.8950 



CM 

Const Slope E? 

146.2816 0.1442 0.6292 

18.10 14.89 
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Figure 2: GM CDS in Blue EDS in Red 
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6 Conclusion 



Approximation of the CGMY Levy measure by hyperexponentials leads to an 
exact Wicncr-Hopf factorization of the approximating process and hence to ex- 
act expressions for the first passage time of the approximating process to a level. 
We then employ these results to obtain closed form formulas for the prices of 
equity default swap contracts that payout on the equity event of a loss of 70% 
of the initial stock value with the receipts being a regular coupon paid till the 
equity event or maturity, whichever is less. 

The methods are illustrated on CGMY processes calibrated to the vanilla 
options market for FORD and GM over the period 25 February 2002 to 25 
February 2005. For the same period we also observe the daily values of the 
credit default swap contracts and these are compared favorably with the option 
imputed equity default swap prices computed by the CGMY approximation 
method proposed here. 
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